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ABSTRACT 


Data previously obtained for the X-33 in the NASA Langley Research Center 20-Inch Mach 6 Air Tunnel have 
been reanalyzed to compare methods for determining boundary layer edge conditions for use in transition correla- 
tions. The experimental results were previously obtained utilizing the phosphor thermography technique to monitor 
the status of the boundary layer downstream of discrete roughness elements via global heat transfer images of the 
X-33 windward surface. A boundary layer transition correlation was previously developed for this data set using 
boundary layer edge conditions calculated using an inviscid/integral boundary layer approach. An algorithm was 
written in the present study to extract boundary layer edge quantities from higher fidelity viscous computational fluid 
dynamic solutions to develop transition correlations that account for viscous effects on vehicles of arbitrary complex- 
ity. The boundary layer transition correlation developed for the X-33 from the viscous solutions are compared to the 
previous boundary layer transition correlations. It is shown that the boundary layer edge conditions calculated using 
an inviscid/integral boundary layer approach are significantly different than those extracted from viscous computa- 
tional fluid dynamic solutions. The present results demonstrate the differences obtained in correlating transition data 
using different computational methods. 

NOMENCLATURE 

h heat transfer coefficient, h=q/(Haw'f^w)^ 


II viscosity (kg/m/s) 

p density (kg/m^) 


(kg/m^/s) 

H enthalpy (J/kg) 

k roughness element height (mm) 

L length of vehicle from nose to end of engine 

module (m) 

M Mach number 

q surface heat transfer rate (W/cm^) 

Re unit Reynolds number ( 1/m) 

Re0 momentum thickness Reynolds number. 


T temperature (K) 

U velocity magnitude (m/s) 

a angle-of-attack (deg) 

5 boundary layer thickness (mm) 

B momentum thickness (mm) 
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Subscripts 

00 freestream static conditions 

aw adiabatic wall conditions 

e local edge conditions 

FR conditions from Fay-Riddell calculation for a 

hemisphere 

w model surface conditions 

tr transition onset 

inc incipient 

ejf effective 

INTRODUCTION 

The ability to accurately predict when boundary 
layer transition will occur on a flight vehicle is impor- 
tant when considering the sizing of a thermal protection 
system (TPS). In recent studies^’^, transition correla- 
tions were generated for a set of wind tunnel data using 
boundary layer edge conditions calculated using invis- 
cid/integral boundary layer techniques. However, it is 
now reasonable to generate transition correlations for 
arbitrarily complex vehicle geometries using boundary 
layer edge conditions extracted from higher fidelity vis- 
cous CFD solutions. This approach was used in a recent 
Mars Science Lander transition study^ and has been 
applied to the X-33 experimental aeroheating database 
in the present study. 
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Berry^ has recently addressed the importance of 
the computational method used when generating tran- 
sition correlations. Different computational methods 
can provide significant differences in the calculated 
edge properties used to form the correlations, as noted 
in Refs, 4 and 5, and the use of a correlation must be 
consistent with the computational method used. 
Therefore, the purpose of the present study is to com- 
pare previous correlations made using an invis- 
cid/approximate boundaiy layer code to correlations 
generated using boundary layer edge quantities 
extracted from viscous CFD solutions. 

The X-33 vehicle, which was intended to be a 
sub-orbital, half-scale representation of a Sin- 
gle-Stage-to-Orbit Reusable Launch Vehicle 

(SSTO-RLV)^’^'^, was selected as the basis for this 
study. Although the program is no longer active, a 
large experimentai/computational aerothermody- 
namic database^" was developed by NASA for the 
X-33, and a considerable investment was made to pro- 
duce computational surface geometries and grids and 
to design and fabricate wind tunnel models. There- 
fore, the X-33 configuration was an ideal choice for 
this study. 

The present study utilized the same discrete 
roughness transition correlation methodologies used 
during investigations into discrete roughness elements 

on Shuttle Orbiter* '*, X-33*', and X-38'^. These 
studies have used the transition parameter Re^M^ (the 
momentum thickness Reynolds number divided by 
the local edge Mach number) along with the boundaiy 
layer thickness, discrete roughness height, and experi- 
mental transition results to correlate the data. 

X-33 VEHICLE GEOMETRY 

The computational results presented in this refer- 
ence are based on the 604B002F configuration of the 
Lockheed-Martin X-33 vehicle, which is commonly 
referred to as the F-Loft, Rev-F configuration. This 
configuration (Fig. 1) is a lifting-body delta planform 
with twin vertical tails, canted fins and body flaps. 


The body length is 19.3 m (63.2 ft.) from the nose to 
the end of the engine module, and the span across the 
canted fins is 23.2 m (76.1 ft.). The canted fins have a 
dihedral of 20-deg and a -8.58-deg incidence angle. 

EXPERIMENTAL RESULTS 

The results presented in Ref. 11, from Test 6770 
in the NASA Langley 20-Inch Mach 6 Air Tunnel, 
were reanalyzed in order to formulate new transition 
correlations based on higher fidelity viscous CFD 
solutions in the present study. Representative flow 
conditions for each of the standard 20-Inch Mach 6 
Air Tunnel operating points at which tests were con- 
ducted have been computed using the GASPROPS^^ 
code and are listed in Table 1. The flow conditions 
listed in Table 1 were used as the free stream flow 
conditions used in the viscous CFD solutions from 
which boundary layer edge conditions were extracted 
in the present study at angles-of-attack of 20-degrees, 
30-degrees, and 40-degrees. 

Global surface heating distributions were 
obtained using the digital optical measurement 
method of two-color, relative-intensity, phosphor 

thermography^ The heating data are presented 
herein in terms of a non-dimensional heat transfer 



Front View 

Figure 1: Dimensioned Sketch of Full-Scale X-33 
F-Loft, Rev-F Configuration. 


Operating 

Condition 

Rejft 


(K) 

P« 

(kg/m^ 

(m/s) 

f^FR 

(kg/m2/s) 

^FR 

(W/cm2) 

«e„ = 9.46xl0^/m 

2.88x10* 

6.00 

62.099 

4.4511x10-2 

946.61 

0.45806 

9.5938 

Re^ = 13.13xl0^/m 

4.00x10* 

6.00 

62.325 

6.1947x10-2 

947.73 

0.54117 

11.403 

Re„= 15.29xl0^/m 

4.66x10* 

6.01 

63.240 

7.2564x10-2 

956.12 

0.59208 

13.004 


Table 1: Wind Tunnel Free Stream Conditions 
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Row Direction 



Top View 



Oblique View 


k St 0.0635 mm 
s 0.1270 mm 
a 0.1005 mm 


Figure 2: IVip Sketch Showing Orientation, 
Width, and Height 



Figure 3: Trip Locations and Fiducial Marks. 

coefficient ratio, h/hfj^, where is the theoretical 

heating computed with the Fay>Riddell^^ method for 
a sphere of the radius of the nose of the test models, 
with a wall temperature of 300 K. 

The cast ceramic aeroheating models were 
0.0 132- scale representations of the full scale vehicle, 
resulting in a 0.254 m (10.0-in.) model length mea- 
sured from the nose to the end of the engine module. 
Discrete roughness elements were produced by appli- 
cation of 0.064 mm to 0.1905 mm (0.0025-in. to 
0.0075-in.) height squares of polyimide film with a 
silicone adhesive. Roughness elements (Fig. 2) fabri- 
cated from the film were applied to the various loca- 
tions of interest on the model and could be easily 
removed without adversely affecting the phosphor 
coating. The roughness elements were placed directly 
over the various fiducial marks both on the center-line 
and on the attachment lines (Fig. 3; see Ref. 1 1 for a 
complete description of test models and discrete trips 
and their placement). 


COMPUTATIONAL METHOD 

Numerical Algorithms 

A boundary layer transition correlation for X-33 

was previously developed^ using the engineering code 
LATCH (Langley Approximate Three-Dimensional 


Convective Heating)^^ to obtain the boundary layer 
edge properties. The LATCH code is an approximate 
three-dimensional heating code based on the axisym- 
metric analog for general three-dimensional boundary 
layers. An integral heating method is used to compute 
the heating rates along three-dimensional inviscid 
streamlines. The inviscid streamlines were supplied 
by the LAURA (Langley Aerothermodynamic 
Upwind Relaxation Algorithm)^’^^ code run in an 
inviscid mode. The LAURA code is a three-dimen- 
sional, finite-volume fluid dynamics solver for 
steady-state flows. Roe-averaging (Ref. 30) with 
Harten’s entropy fix (Ref. 31) and Yee’s Symmetric 
Total Variation Diminishing limiter (Ref. 32) is used 
for inviscid fluxes, and a second-order scheme is 
employed for viscous fluxes. Further details regarding 
the LATCH computations can be found in Ref. 10. 

For the current study, laminar, perfect-gas air 
aeroheating computations were performed at the wind 
tunnel test conditions using the thin-layer 
Navier-Stokes LAURA option. Free stream conditions 
for the computations were set to the nominal operat- 
ing conditions of the NASA LaRC 20-Inch Mach 6 

Air Tunnel, which are listed in Table 1. For these com- 
putations, a uniform, ambient 300 K wall temperature 
boundary condition was imposed. The use of a con- 
stant wall temperature was valid because the experi- 
mental data were reported in terms of the 
non-dimensional ratio, /i//i/r/^, which is approximately 
constant with wall temperature. The computational 
grid used in Ref. 12 was utilized in the present study 
(see Fig. 4). This grid, which was originally used in 



Figure 4: Bow shock aligned computational grid 
for vehicle at a= 30-degrees. 
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the GASP^^ CFD solver, contained approximately 
1.73 million grid points. The grid was decomposed 
into 16 computational blocks and was adapted to align 
the grid with the bow shock and cluster grid points in 
the boundary layer with a wall cell Reynolds number 
of approximately 3. 

A comparison of windward centerline heating 
distributions between computations using LAURA 
and experiment at a = 20-degrees, 30-degrees, and 
40-degrees is shown in Figure 5. The experimental 
data, which are taken from Ref. 1 1, show an increase 



xA. 



0.00 0.20 0.40 0.60 0.80 1.00 

XA. 



0.00 0.20 0.40 0.60 0.80 1.00 

x/L 

(c) a s 4O-d0gr»»s 

Figure 5: Comparison of Measured (see Ref. 11) 
and Predicted Values of Centerline Heating. 


in heating near the aft region of the model as the free 
stream Reynolds number was increased. This is due to 
the onset of natural transition. The laminar computa- 
tional predictions generally agree within ±10% of the 
laminar experimental data at all values of Re^ (there 
is no Reynolds number effect on non-dimensional 
laminar heating). 

Boundary Layer Edge Definition 

Engineering relations^*^ are used within LATCH 
to obtain both the momentum thickness (6) and the 
boundary layer thickness (8)^ where 6 is calculated 
based on local properties and 8 is determined based 
on the shape factor relation of 8/6 = 5.5. The bound- 
ary layer edge properties are obtained by interpolating 
in the inviscid flowfield a distance equal to the bound- 
ary layer thickness away from the wall. To accomplish 
this, an initial assumption is made for the boundary 
layer edge properties (usually the wall values), and 
the boundary layer thickness is computed. Then the 
edge properties are re-computed based on this new 
location within the flowfield and the solution is iter- 
ated until the re-computed boundary layer thickness is 
equal to the assumed value. This process usually takes 
two or three iterations to converge. The use of edge 
properties determined in this manner approximately 
accounts for the effect of variable entropy at the 
boundary layer edge. Recent experience with moder- 
ately blunt bodies like the Orbiter, however, has 
shown that based on the inviscid solutions available, 
iterating away from the surface has a minimal affect 
on the convective heating solutions (the primary moti- 
vation for the code), even though it provides for more 
exact edge conditions. By assuming the boundary 
layer edge properties to be equal to the inviscid sur- 
face conditions, which corresponds to a constant 
entropy condition, solutions are obtained much more 
quickly. Since the correlations previously formed for 
the X-33 were based on a constant entropy assump- 
tion, the boundary layer edge conditions from the 
inviscid/integral results presented herein are actually 
the inviscid wall conditions. 

In order to extract boundary layer edge condi- 
tions from the LAURA solutions, a post-processing 
algorithm was written. This algorithm computes sur- 
face normals at each computational point on the sur- 
face of the model and iterates along the normal until 
the boundary layer edge is found, which is defined as 
a specified value of the ratio (discussed later). 

Once the edge is found, flow properties are interpo- 
lated at discrete points along the surface normal in 
order to extract the boundary layer profile. The 
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momentum thickness can then be calculated by inte- 
grating over the boundary layer. 

The criteria for which the edge of the boundary 
layer is defined can change throughout the speed 
regime from subsonic through hypersonic. For many 
flows, the boundary layer edge is defined where the 
velocity is equal to 99% of the free stream velocity. 
However, this assumption is not always valid for 
hypersonic flows, where viscous effects can be much 
more significant. The criterion most often used to 
define the boundary layer edge in hypersonic flows is 
the location along the wall normal at which the local 
stagnation enthalpy (//^) is equal to 99.5% of the free 
stream stagnation enthalpy. This definition can be 
used because, for an inviscid, adiabatic steady flow 
with no body forces, the total enthalpy is constant 
along a given streamline. Thus, outside of the bound- 
ary layer, the stagnation enthalpy should be constant. 

BOUNDARY LAYER EDGE PROPER^ 
TIES COMPARISON 

Sensitivity to Boundary Laver Edge Definition 

The sensitivity of the boundary layer thickness, 
momentum thickness Reynolds number, and edge 
Mach number with respect to the definition of the 
fraction of the local stagnation enthalpy to the free 
stream stagnation enthalpy are shown in Figures 6, 7, 
and 8, respectively. Data are shown along the wind- 
ward centerline at angles-of-attack of 20-degrees, 
30-degrees, and 40-degrees and a fixed free stream 
Reynolds number of 13.13xl0^/m. Total enthalpy 
ratios of Hq/Hq oo = 0.950, 0.990, and 0.995 are pre- 
sented from the post-processed LAURA solutions. 
The Hq/Hq 00 = 0.950 boundary layer height, momen- 
tum thickness Reynolds number, and edge Mach num- 
ber calculations were approximately 20%, 30%, and 
10% lower than the Hq/Hq oo = 0.990 calculations, 
respectively. The difference between the Hq/Hq = 
0.990 and 0.995 calculations is below 5% for the three 
boundary layer properties considered. This suggests 
that defining the boundary layer edge as the location 
where the stagnation enthalpy ratio is Hq/Hq ^o = 
0.995 is appropriate. 

Comparison of Method of Calculation 

Comparisons of boundary layer edge parameters 
along the windward centerline between LATCH solu- 
tions and the post-processed LAURA solutions are 
also shown in Figures 6-8. The windward centerline 
boundary layer height, 6, is presented in Figure 6. For 
all angles-of-attack, <5 at the nose is comparable 


between the LATCH and post-processed LAURA 
solutions. However, as the running length increases, 
the boundary layer height extracted from LAURA 
increased at a greater rate than the LATCH solutions. 
This trend decreased as angle-of-attack was increased. 
Therefore, it appears that the LATCH solution was 
under-predicting the boundary layer height, especially 
at a = 20-degrees. This greater difference at 
20-degrees may be due to the thickening of the 



Figure 6: Boundary Layer Thickness as Calcu- 
lated from LAURA at Varying Stagnation 
Enthalpy Ratios and from LATCH at Reoo = 

13.13xl0*/m. 
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boundary layer along the centerline due to inflow 

from the attachment lines^^ 

For Rc 0 (Figure 7), the LAURA and LATCH 
solutions compared favorably up to approximately x/L 
= 0.25, but the solutions diverged downstream of this 
point, with the LAURA solution predicting a higher 
value of Re 0 . Once again, the difference between the 
two solutions became greater as angle-of -attack was 
decreased. 



(a) a * 20-degr999 




Figure 7: Momentum Thickness Reynolds Num- 
ber as Calculated from LAURA at Varying Stagna- 
tion Enthalpy Ratios and from LATCH at = 

13.13x10*^/111. 


The divergence of the computed values of edge 
Mach number (Figure 8), however, was in the oppo- 
site direction. The post-processed LAURA solutions 
predict a value of which was lower than that of the 
LATCH solution. The edge Mach number calculated 
by post-processing the LAURA solution was approxi- 
mately 20% lower than that calculated by LATCH. 

When comparing methods of calculating bound- 
ary layer edge conditions using inviscid/engineering 
methods. Berry ^ showed that when the boundary layer 




0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 

X/L 

(c) a s 4C-d«gr»€s 


Figure 8: Edge Mach Number as Calculated from 
LAURA at Varying Stagnation Enthalpy Ratios 

and from LATCH at Re« = 13.13xl0^/m. 
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transition parameter Re^M^ was examined, the meth- 
ods fortuitously showed good agreement, due to the 
cancellation of differences in Rcq and M^. This 
parameter is shown along the windward centerline at 
a - 20-degrees, 30-degrees, and 40-degrees for the 
LATCH and post-processed LAURA {Hq/Ho oc - 
0,995) solutions in Figure 9. The two solutions com- 
pare well upstream of approximately x/L = 0.25. 





Figure 9: Re^Mg as Calculated from LAURA at 
H^Ho co = 0.995 and from LATCH at Reoo = 
ISJaxlO^^/m. 


Downstream of x/L = 0.25, however, the transition 
parameter calculated by LATCH was lower than that 
of the post-processed LAURA solution. This differ- 
ence once again became more pronounced as 
angle-of -attack was decreased. 

BOUNDARY LAYER TRANSITION 
CORRELATION 

The results presented in Ref. 1 1, from Test 6770 
in the NASA Langley 20-Inch Mach 6 Air Tunnel, 
were used to generate a new, higher fidelity transition 
correlation based on boundary layer edge conditions 
extracted from the viscous CFD solver LAURA. 
Approximately 100 data points were used from this 
test in the current study. Data points from along the 
centerline and along the attachment lines were 
included. Laminar computations were performed 
using LAURA for a complete wind tunnel test range 
of angles-of-attack (a = 20-deg, 30-deg, and 40-deg) 
and a selection of Reynolds numbers (Rcao = 

9.46xl0^/m, 13.13xl0^/m and 15.29xl0^/m) in order 
to determine boundary layer edge properties (Re^, M^, 
and 5), which were used to formulate transition corre- 
lations from the experimental data. 

For each discrete roughness element in each 
experimental run which was reanalyzed, the state of 
the boundary layer downstream of the element was 
determined through visual inspection of the surface 
heating images and classified as either: laminar; tran- 
sition at some distance downstream of the discrete 
roughness element; or transition effective (i.e. imme- 
diately downstream oQ from the discrete roughness 
element. Examples of experimental data which fit 
each of these classifications are shown in Fig. 10. 

The results of X-33 transition correlations gener- 
ated using LATCH were presented in Ref. 2. Two sim- 
ple curves were suggested as laminar and turbulent 
boundaries (shown in Fig 1 1 as dashed curves). The 
first curve (incipient) was defined as the boundary 
between laminar flow downstream of the discrete 
roughness element and transitional flow downstream 
of the element. The second curve (effective) was 
defined as the boundary between transitional flow 
downstream of the element and where the transitional 
wedge was attached to the roughness element. 

Data points taken from the experimental study 
where the boundary layer edge conditions were 
instead extracted from LAURA solutions are also 
shown in Fig. 1 1. The curves suggested in Ref. 2 gen- 
erated using the LATCH edge conditions were not 
applicable when using the LAURA edge conditions, 
as evident from the discussion on boundary layer edge 
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Figure 10; Classification of Discrete Roughness 
Effect on Boundary Layer. 


properties earlier in this paper. This does not mean 
that the LATCH curves are incorrect. However, the 
present results demonstrate the differences obtained in 
correlating transition data using different computa- 
tional methods. Based on the results of this study, two 
simple curves similar to those of Ref. 2, but for the 
LAURA boundary layer edge conditions, are sug- 
gested in Figure 11. 

CONCLUSIONS 

Data previously obtained for the X-33 in the 
NASA Langley Research Center 20-Inch Mach 6 Air 
Tunnel have been reanalyzed to form new boundary 
layer transition correlations using higher fidelity vis- 
cous computational fluid dynamic solutions. An algo- 
rithm was written in the present study to extract 
boundary layer edge quantities from viscous compu- 
tational fluid dynamic solutions to develop transition 
correlations that account for viscous effects on vehi- 
cles of arbitrary complexity. 



0.00 0.20 0.40 0.60 0.80 1.00 1.20 

Figure 11: LATCH Transition Correlation Curves 

with Data Points from LAURA Solutions. 

It was found that defining the boundary layer 
edge as the location where the ratio of local stagnation 
enthalpy to free stream stagnation enthalpy is equal to 
0.995 was appropriate. When the LATCH and 
LAURA boundary layer edge quantities were com- 
pared, the boundary layer height and momentum 
thickness Reynolds number from LAURA were sig- 
nificantly higher than the inviscid/integral values; 
whereas the viscous edge Mach number was lower 
than the computed inviscid/integral value. These dif- 
ferences were found to increase as angle-of-attack 
was decreased. 

When the transition correlation Re^M^ vs. k/5 
was revisited using the boundary layer edge quantities 
extracted from LAURA solutions corresponding to 
wind tunnel conditions, a different set of transition 
curves were required. This does not imply that the 
correlation generated using the LATCH boundary 
layer edge quantities was incorrect, only that when 
one method is used to obtain boundary layer edge 
conditions, the corresponding set of curves must be 
used. 
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